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The Magnetized Electron Gas in terms of 
Hurwitz Zeta Functions 



We obtain explicit expressions for thermodynamic quantities of a rela- 
tivistic degenerate free electron gas in a magnetic field in terms of Hurwitz 
Zeta functions. The formulation allows for systematic expansion in all 
regimes. Three energy scales appear naturally in the degenerate relativis- 
tic gas: the Fermi energy E F , the temperature T and an energy related 
to the magnetic field or Landau level spacing, eB/E F . We study the cold 
and warm scenarios, T <C eB/E F and eB/E F <C T, respectively. We re- 
produce the oscillations of the magnetization as a function of the field in 
the cold regime and the dilution of them in the warm regime. 
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1 Introduction 



The study of the thermodynamic properties of degenerate relativistic electron 
gases in strong magnetic fields, as those found in compact astrophysical objects, 
was started long ago Although in that work the problem was stated in 
quite general terms, that is, arbitrary temperature and relativistic electrons 
described by Dirac's equation, approximate analytical results were obtained only 
in the limits of low temperatures and non-relativistic electrons. While the non- 
relativistic approximation has some validity in white dwarfs, it is inappropriate 
for the electron gas existing inside neutron stars. In fact, in a typical neutron 
star the electron density is of the order of 10~ 2 fm~ 3 , which implies a Fermi 
kinetic energy T F = E F — m of the order of 100 MeV. The temperature, on the 
other hand, is at most a few MeV, which makes the ratio T/T F a small parameter 
and the gas quite degenerate. The relativistic case was partially considered 
later on in Refs. |2[^| with regard to the problem of self-magnetization, and in 
Ref. |Q in connection with the magnetic susceptibility. In all these works only 
the dominant contributions relevant for each case were kept in the analysis. 

More recently, a renewed interest in the relativistic magnetized electron gas, 
both at zero and finite temperature, has appeared from the point of view of 
quantum field theory [||-^). In these studies most of the emphasis is put on 
the formal aspects of the problem. For instance, in Ref. [Q an exact analytic 
expression is obtained for the effective action (also called the grand potential 
in statistical mechanics), which corrects an incomplete result given earlier || 
that had missed the de Haas - van Alphen oscillations that are present at low 
temperatures, exactly as in the non-relativistic case. In Ref. || the quantum 
field theory results of Ref. [Q are recast, for the case of low temperatures, in 
terms of sums over filled Landau levels, which is actually the starting point of the 
quantum statistical mechanical calculation developed in the earlier works Jl|-Q . 

In this work we revisit the usual elementary statistical mechanical approach, 
originally carried out by Landau fl(],[ll| (see also Ref. JlJ]) for a non-relativistic 
gas, to study a highly degenerate relativistic free electron gas, in the presence 
of a uniform magnetic field. Starting from a closed analytical expression for the 
density of states for this system, we derive explicit, simple analytical expressions 
for the various quantities of thermodynamic interest, such as the density and 
magnetization. The object of central interest will be the grand potential, from 
which all relevant thermodynamic quantities can be computed. 

It is well known that for a relativistic electron gas the Coulomb interactions 
among electrons and between the electrons and any background of positive 
charge that should exist in a neutral system are small (order a) corrections 
to the kinetic energy, so we will neglect them. We will therefore work consis- 
tently to zeroth order in QED corrections. In particular, this means that we will 
take the electron gyromagnetic ratio exactly equal to two, g — 2. Additionally, 
since we shall be solely concerned with the degenerate regime, we neglect the 
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contribution of positrons altogether. 



We will show that the grand potential at T — has the integral representation 

y f (E F - m 2 )/2eB n J\ 

Q (E Fl B) = - JL (2eB)^ / 7=T== d «> W 

^ Jo \/m 2 + 2eB q 

where m and e are the electron mass and the fundamental charge, respectively 
(we use natural units, so that % = c = ks = 1)- Consequently the electron 
density at T = is 

The function Tt -1/2(9) appearing in both (0) and (§), is a combination of Hur- 

witz zeta functions, H- 1/2(9) = C(~5j {9})~C(~ 51 9+1) — 5"\/^' anc ^ pl & y s a per- 
vading role in our study of the magnetized free electron gas. Since the pressure 
is given by P = —Q/V, the expressions (|l|) and @) furnish a parametric repre- 
sentation of the equation of state at T = 0, Pq — Po(E F , B),uq = no(E F , B). 

The expressions given above correspond to zero temperature, but, as we show 
later in the paper, for a system of non-interacting fermions the full finite tem- 
perature grand potential can be obtained from its zero temperature counterpart. 
For instance, in the degenerate regime finite temperature corrections are usually 
obtained systematically from the ground state quantities in the form of a Taylor 
series in powers of T//i, where /i is the chemical potential. Take, for example, 
fl(T,fi). As is well known, the leading correction to the zero-temperature re- 
sult generally goes as (T/E F ) 2 . This will actually be the case, provided the 
grand potential at T — 0, Qq(E f ), does not vary greatly when E F is changed 
by an amount of order T. However, as we will see, in our case Qq(E f ) has 
contributions that oscillate rapidly with E F due to the filling of discrete Lan- 
dau levels. Therefore the approximation above is not valid, except for very low 
temperatures (T eB/E F , as we will see later). We will compute with our 
formalism the correct finite temperature behavior of Q(T, fi) and reproduce the 
previously known result that the oscillations are smoothed out as the tempera- 
ture increases. 

The novelty of our approach is the use of the Hurwitz Zeta function to deal 
with the discrete sums over Landau levels that accommodate the electrons in 
the presence of a uniform magnetic field. The traditional approach [^2| makes 
use of the so-called Poisson summation formula which is of limited use if one 
needs to expand the resulting expressions in different regimes, such as for small 
fields. The Hurwitz zeta function approach is much more powerful in this case, 
as it leads to closed analytical expressions, which can eventually be evaluated 
numerically with ease, and is suitable for expansion in any desired regime. Al- 
though this is not the first appearance of the Hurwitz zeta function in connection 
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with the system being studied MM, this seems to be the first time that its an- 
alytic properties are fully put into use to unravel the thermodynamics of the 
relativistic magnetized free electron gas. 

The Hurwitz zeta function £(z, q) is defined as the analytic extension to the 
whole complex z plane of the series 

valid for Hez > 1 and q 7^ 0, — 1, — 2, The resulting function is analytic 
everywhere except at z = 1, where it has a simple pole with unit residue (see 
Appendix A for details). 

In section 2 we formulate the density of states in its most general form in terms 
of Hurwitz zeta functions. In section 3 we classify the physical regimes we 
want to study. The core of the thermodynamics is presented in Section 4 and 
particular studies of the magnetization are shown in Section 5. Details of the 
calculations and mathematical formulas are given in the appendices. 



2 The Density of States in terms of Hurwitz 
Functions 

The stationary states of a Dirac electron that moves in a uniform magnetic field 
B (which we take to point along the z-direction), are specified in terms of four 
quantum numbers. In the gauge where A(x) = x B y, these quantum numbers 
are p z (momentum component along B), p y (momentum component along A), 
a non-negative integer n = 0, 1, 2, . . . that specifies the Landau level, and an 
integer A = ±1 which denotes the spin parallel or antiparallel to B. The energy 
levels are independent of the quantum number p y and are given by^ [|l3| 

E(p z ,j) = ^/pl+m 2 + (2eB)j, (4) 

where j is a nonnegative integer defined by j — n + (A + l)/2. The energy 
levels (|J) are highly degenerate, due to their independence on p y and their 
dependence on n and A only through the single combination j. The number of 
states g(p Zl j)dp z for given j and momentum along B between p z and p z + dp z , 
in a system of electrons confined to a finite cubic box of volume V, is given by 

eB 

g(Pz,j)dp z = gj^Vdp z , (5) 

with gj = 1 for j = and g* — 2 for j ^ 0. Physically, this degeneracy embodies 
the fact that the levels with j = accommodate only electrons with spin point- 
ing down, whereas all the others can have electrons spinning in either direction. 

3 For g = 2, i.e., neglecting QED corrections to the electron's magnetic moment. 
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As an immediate consequence, the levels with j = alone are responsible for 
the paramagnetic contribution to the magnetization. See [Eij for details. The 
density of states in energy space is thus given by 

JE 

g(E)dE = 2j29(\p*\,j)d Pz , (6) 

3=0 

where the factor of 2 in front of the sum takes into account both possible signs 
of p z , and the integer j E corresponds to the highest Landau level j that starts 
at an energy less than E. From Eq. (^) it is clear that j E is the integer part of 
the quantity 

E 2 -m 2 

* = -£b-> (7) 

denoted usually as j E = [q E \ . The degeneracy of the levels |l3| can be expressed 
explicitly by 



£ y/E 2 -m 2 - 2eB ■ j 



j 

L«eJ 



47r f-^ V1e-J 

Considering now the degeneracy gj in Eq. (||) and formula ( |6l| ) with z = 1/2 
we find 



L«eJ 

E 



9j- 



]=0 V9e-J 



C(i &- LfeJ)-C(ife + i) 



2\/9e 



(9) 



where the r.h.s. of this equation defines the following family of functions that 
appear regularly in the thermodynamic expressions for this system: 

H z (q) = C(z, {q}) - C(z, q+l)- \f\ (10) 
and where {q} — q — [q\ is the fractional part of q. 

In what follows, we will denote energies in units of the electron mass, e = E/m, 
and the magnetic force eB in terms of m 2 : 

, 2eB 



(11) 
iii- 

The variable b measures the magnetic field strength in units of the natural 
strength B = ±m 2 c 4 /ehc w 2.2 x 10 13 Gauss §J|. 

With this notation and the use of Eq. (||) , the density of states for energy E has 
a simple form: 

9(E) = v£&* s H 1/2 (^p) . (12) 
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The total number of states up to a given energy E is the integral over the density 
of states: 

G{E) = f g{E')dE'. (13) 

Using Eq. ( |62| ) of Appendix A, it is easy to show that W 1 / 2 ('z) = 2H'_ 1 y 2 (<7) 
(where the prime denotes differentiation with respect to q), so that we can do 
this integral in a closed form, obtaining: 

G(E) = V £** «_ 1/a (^) . (14) 

Clearly, G(E F )/V is the electron density in the system, where E F denotes the 
Fermi energy. We will consistently use dimensionless variables in the expres- 
sions: 



e F = — , Pw = y/e F -l. (15) 

771 

While e F is the Fermi energy, p F is the Fermi momentum of the j — Landau 
level, in units of m. The expression p F /b gives a measure of the occupation of 
the Landau levels in the ground state: its integer part, [p F /b\, equals the last 
Landau level that is occupied, and its fractional part {p F /b}, or more precisely 
1 — {p F /b}, gives a measure of the "distance" to the next empty one. 



3 Limiting cases and regimes 

In the non-relativistic case, the energy scale associated with the magnetic field 
is the spacing between Landau levels, eB/m, which is uniform and equal to ojb, 
the cyclotron frequency. For typical values of magnetic field strengths found in 
pulsars, say B = 10 12 Gauss, one finds lub ~ 10 keV. 

As is well known, a fermion gas becomes degenerate when the occupation num- 
ber falls abruptly for states with energy above a given value E Fl the Fermi 
energy. This happens if the temperature is much smaller than the kinetic range 
T F = E F — m {i.e. the energy range between the lowest orbital and the Fermi 
energy). For electrons in neutron stars, typical values might be T < 1 MeV 
and T F ~ 10 2 MeV, so that the gas is not only degenerate but also highly rela- 
tivistic. In the relativistic regime however, the level spacing is not uniform, but 
decreases with energy, so that close to the Fermi surface it is eB/E F . 

While the Fermi energy is invariably larger than the temperature in a degenerate 
gas, the number of occupied Landau levels may or may not be large. This 
number is given by the integer part of (-E 2 — m 2 ) / (2eB) , in view of Eq. (g) and 
its preceding paragraph, and so it is large if 2eB/E F <C T F . This is indeed the 
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case for typical magnetic fields in neutron stars (B < 10 14 G) where eB/E F < 10 
keV. 

We can still classify the gas in a uniform magnetic field as cold of warm, de- 
pending on whether the temperature is smaller or larger than the energy spacing 
between Landau levels close to the Fermi surface, the latter being eB/E F . Con- 
sequently, in a cold degenerate gas, the following condition is satisfied: 

eB 

T <C — and T <C T P . (16) 
E F 

If, in addition, the magnetic field is small so that there is a large number of 
occupied Landau levels, the above hierarchy becomes: 

T « ^ « T F . (17) 
E F 

On the other hand, in the warm regime, the temperature is comparable or larger 
than the Landau spacing: 

^ < T « T F . (18) 
E F 

These three cases are the regimes of interest here. 



4 The Grand Potential and Hurwitz Functions 

The general expression for the grand potential for an ideal Fermi gas in a uniform 
background magnetic field B is given by 

Sl(T, B, V, M) = -T E ln i 1 + e-^-^ T ) , (19) 

a 

where /x is the chemical potential and the sum is over all 1-particle orbitals a. 
The B dependence is implicit in the orbital energies E a and the density of states. 
Thermodynamic quantities can be obtained directly from f2. In particular, M., 
the magnetization per unit volume is obtained as 

In terms of the density of states g(E) the grand potential (per unit volume) is 
oj = — = -T f^ln(l + e -(E-n)/T\ dE ^ (21) 

This expression, together with the explicit form ( |l^ ) for g(E), can be used to 
study the thermodynamics at arbitrary temperature T and magnetic field B. 
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The T = Limit 



As we show in Appendix B, for a system of non-interacting fermions any finite 
temperature quantity can be obtained from its zero-temperature counterpart. 
Therefore, in the rest of this section we will just concentrate on the grand 
potential density at zero temperature, u>(T = 0) = u)q, keeping in mind the 
standard definition for the Fermi energy E F = fi(T = 0). The T = limit is 
simply formulated from 



lim 



Tin 1 



-(E-ix)/T\ _ 



E 




for E < E F 
for E> E v 



so that 



pE F 

uj (E f ,B) = / 

J m 



9(E) 
V 



E dE- E F 



N 

V 



(22) 



(23) 



The integral in Eq. (]25|), which is the ground state energy density (= Uq), can be 
done using the explicit representation for the density of states given in Eq. ( |l2| ) . 
After standard manipulation and integration by parts we get: 



U 



E F 



9(E) 
V 



EdE = 



m 
4^2 



2 6 3 / 2 e F H-x/2 



^ r %/b ^M dq 



(24) 



Notice that the first term above is exactly E F N/V, according to Eq. (|14[), so it 
cancels with the second term in Eq. (p3|), leaving a compact form for the grand 
potential density u>q = Qq/V: 



w (e F ,b) 



:b 5 ' 2 



Pi/ b H- 1/2 (q) 



dq. 



where 



H_ 1/2 (g) = C(-|,W)-C( 
Formally, this is the main result of this work 



hi 



(25) 



(26) 



Expression (|2^) allows us to recover the well known non-relativistic result |L2| 
in a straightforward fashion: in the non-relativistic limit p F -C 1, so that the 
square root in the integrand of (25) can simply be replaced by unity. Using then 
the relation 7i-i/2(q) 



1^-3/2 (<z) w e find 



(n.r.) 



(e F ,b) 



6tt 2 



b 5 / 2 H_ 3/2 (p 2 F /b) 



where 



w_ 3 /2(g) 



C(-1,W)-C(-1 



,9+1) 



I 3/2 



(27) 



(28) 



8 



The result (|27|) is exactly the known non-relativistic result, after identifying 
pi = 2/i„. r ./m, where fi n . r . is the standard (non-relativistic) chemical potential. 
In particular, notice that Ti._ 3 ^ 2 (q) splits naturally into a monotonic term and 

an oscillatory term, 7^—3/2 (?) = T^L^il) +^-3/2(9)1 as: 

W^(9) = -CK,<Z + l)-^ 3/2 , (29) 
H%{q) = C(-|,W)- (30) 

The oscillating character of 7i_l c ) 2 (Q) arises from the fractional part {q} in the 
definition above. This is the term that accounts for the de Haas - van Alphen 
oscillations in the non-relativistic case, and leads to Landau's result at finite 
temperature, as it is shown in section ^. 

Likewise, the relativistic grand potential density (^5| ) separates into two terms, 

(moti) , (osc) n 

ojq = uiq + luq , where 

^, b) _ ^f* <ti££M «, (31) 

Here we have kept the notation consistent, however one must be aware that 
loq ^ , the integral of an oscillating function, is not purely oscillatory, as we will 
see. The integral expressions above are exact and finite for any values of the 
magnetic field and Fermi momentum. 



Although the integrals in Eqs. ( |3l| ) and (32) cannot be solved in a closed form, 
one can nevertheless use the analytic properties of the Hurwitz zeta function 
to obtain an expansion of too for small 6, or similarly for large "p\jb (i.e. large 
number of occupied Landau levels). Let us first consider WQ mon \ A naive bi- 
nomial expansion of the square root in the integrand will not work because of 
the diverging behavior of the integral as p\jb —> 00. A procedure to obtain the 
correct expansion for small b consists in first extracting out a sufficient number 
of leading terms in the asymptotic expansion of the zeta function for large q, 
such that the remainder is integrable in the limit b — ► [three terms will do 
for £(— 1/2, (?)], and then integrating the subtracted terms explicitly. Based on 
Eq. ( fn] ) of Appendix A we write 

C(*. q) = ^tT 2+1 + + ^q-'- 1 + g), (33) 

implicitly defining the subtracted zeta function (^3(2;, q). Now, the integrals cor- 
responding to the exhibited terms in (^3|) along with the term yfq/2 in ( |3l| ) can 
be done exactly for finite &, leaving only an integral containing the subtracted 
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function fa(—l/2,q): 



(mon) 



i 



1. 



4tt 



G 



Vb 



,V2 



1 

2 

1 r 
3 

[cosh _1 (£ F ) -_p F e F ] 



(b + p 2 F ) ' e ¥ 
(b + plf 2 * 



1 

6 3/2 



6 5/2 



pg - /b Cs(-|,g + 1) 



(34) 



The expression above is still exact, but now the integral containing £3(2,9) 
can be expanded in powers up to order b without entering into trouble, because 
£;(— 1/2, q) vanishes as q~ 5 ^ 2 for q — > 00 (see Eq. ( |7^ ) of Appendix A). It suffices 
to take the upper limit to infinity and expand the denominator: 



VTT69 



-('?c 



1 

GO 



(35) 



1260 



+ 0(fo 3 / 2 ) 



This integral thus loses the dependence on the chemical potential and becomes 
a spurious contribution that cancels in luq when we put all pieces together (see 
below). A word of caution is due: to go to higher orders in b one cannot 
continue this expansion in the same fashion. Instead, one should follow the 
same procedure from the beginning, but extracting more explicit terms from 
the asymptotic expansion (|33"|), integrating them directly, and then expanding 
the n-residual integral containing ( n (—l/2,q) up to the last convergent term, 
that goes as b 2n ~ 5 . 

Now, let us turn to the oscillatory contribution u;q OSC ^ , which can be written in 
the following form (see Appendix C): 



(osc) 



4 

m 


b 2 


~4^ 








+/ 


a 


Jo 





{pl/b} 



C(-i?) 



y/q + l/b+\p*/b\ 



dq 



(36) 



whose expansion for small b is straightforward (however tedious): 



(osc) b-+0 




m 4 f 2 
4^2 I 3 



6 5/2 



~cH,{rf/6})-c(-§) 



(37) 



10 



15 



-b 7 ' 2 



rC(-|.W/6»-C(-f) 



0{b*' 2 ) . 



Finally we can put all the terms together and find the expansion of the grand 
potential density wo for small magnetic field (see Appendix D): 



™> 1 , -1/ x 1 3 1 
UJ ° = ^4^2 12 ( £ f) + g£ F - Pf 

b 2 2 b 5 / 2 

+ ±- C0S h-\e F ) + l b —C(-l,{p 2 F /b}) 



2 b 7 ' 2 
15 e| Cl 



{jg/6}) + 0(& 4 ) 



(38) 



We note that expansion ( pq ) is in agreement, apart from an overall sign, with 
the small field result for the effective action obtained in Ref. || (Eq. (74)), which 
in our notation reads 



AS d 



=31 



+ ■ 



bra 1 

T=0;/i^0 ~~ "g^2 

6 2 m 4 ^ (b/el) l +V 2 



e F p F 



2ir 2 



E 

1=0 



(2Z + 1)(2Z- 



{ 

3){« 



ln(e F + p F )| 



3 
2 ' 



3 „2 
2 ' "f 



/&)} (39) 



Although the non-analytic oscillatory terms are immediate to compare, it is 
quite nontrivial to obtain the analytic terms from (|39|), since each term of the 
type b l+1 / 2 ((~l - l,p 2 /b) in the sum in ( |39|) c ontributes to all orders in b as 
b — > 0, in view of the asymptotic expansion (72|) for the Hurwitz zeta function. 
Using the first three terms of this asymptotic expansion, one can perform each 
of the resulting infinite sums in a closed way to indeed obtain the terms shown 
in ([38]), with a vanishing coefficient for the term linear in b. In retrospect, we 
realize that the result (39) can be formally obtained directly from the closed 
expression ( p5| ) by doing a binomial expansion of the square root and then 
integrating term by term. However, the latter expansion is valid over the whole 
integration range only if p 2 < 1, i.e. basically in the non-relativistic limit. 



From a numerical point of view it is clearly advantageous to use expansion 
instead of (39), since then one can compute to any desired precision at small b 
by keeping only a finite number of terms in the expansion. 

All thermodynamic quantities at T = can be obtained from loq as derivatives. 
One should notice that, since C( z , 0) = ((z, 1) for z < 0, the function £(z, {q}) is 
continuous in q for all z < and, in view of Eq. (|62|), has continuous derivative 
for all z < —1. For example, it is straightforward to derive the expansion for 
the density at T = using the thermodynamic identity 



n 



1 du>o 
m de F 



(40) 
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We thus find the the density at T — 0, consistent with Eq. (|l4]), and the corre- 
sponding expansion for small 6: 

n (e Fl b) = £ b 3 2 H J 



2tt 2 

m 3 f 2 

2^2 



^Pl + bV2 C(-h{pl/b}) + ^ FF + 0(b%j (41) 

In the b = limit we recover the free electron gas result, no = m 3 p F /3ir 2 . As 
seen in Eq. (pi]), the leading correction for finite b and T — is a term that 
oscillates with p F /b, corresponding to the de Haas- van Alphen effect for the 
magnetization in metals. 

Fig. 1 shows no vs. e F for two values of the field. The step-like behavior is due 
to the filling of consecutive Landau levels as e F increases. The density of states 
in a given Landau level [see Eq. (||)] goes as bdp z ~ bdE/p Zl for p z starting 
from zero. The higher density at the bottom of each level causes the step in 
no- On the other hand, as b decreases the steps gradually disappear until the 
smooth 6 = limit is reached. Alternatively, Fig. 2 shows no as a function of the 
field 6, for fixed Fermi energy. Here no oscillates as 6 changes. Imagine we start 
from a large value of 6. As 6 decreases, the Landau levels move down in energy, 
crossing one by one the threshold defined by e F . Since the levels are denser at 
the bottom, n first grows as each level becomes accessible, but then decreases 
again, because of the 6 factor in the density of states. The oscillation amplitude 
becomes smaller and smaller as the Landau levels get closer and closer, until we 
reach the 6 = limit. 

We note that the expansion in Eq. (|4|) is actually a large p F /b expansion, 
appropriate when a large number of Landau levels are occupied at T = 0. So 
it will remain valid even for large values of the magnetic field, say 6 ~ 100, 
provided the ratio p F /b stays much larger than one. This will be the case for 
electron densities in typical neutron stars, for which p\ ~ 10 4 . 

We can also obtain the expansion for uo, the energy density at T = — whose 



exact integral form is given in Eq. (24) — , using the expansions in Eqs. ( p8[ ) and 
( ^l| ) and the thermodynamic relation uo = loq + m e F no'. 

m A \ 3 1 1 , -1/ \ 

+26 3 / 2 £f C (-1, {Pl/b}) + - J£ cosh- 1 ^) 

2 6 5 / 2 2 6 7 / 2 

-l—ci-IApl/^-^cHAm) 

+ 0(6 4 )...). (42) 
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Unlike Eq. ( pl| ) for the particle density, the expansions for the energy density 
given above and for the grand potential given in Eq. d38| ) are valid for small 
b only. More useful in the case of neutron star conditions is an expansion for 
large p F /b, that is, many Landau levels occupied, but regardless of b. In that 
case, for instance, the small b expansion of the denominator in Eq. (|3^) is not 
valid. The correct expansion of ujq for large p F /b is the following, where some 
p F -independent integrals are left to be done numerically: 



(mon) 



m 4 








1 

~ 2 


> 


1 


\b 



e F + y/b + p F 
1 + Vb 



{b + pl) 
{b + p 2 F ) 



1/2 



3/2 



fel/2 
5 3/2 



+ -b [cosh 1 (e F ) -p F e F ] 
J VTTb^ 



1 fr 4 
3840 si 



(43) 



(osc) -— J -— C(-|,{^M) + 4^rC(-I,{^}) 



4tt 2 3 e F ^ 2 ' l ^ F/ J; 15 e 



+ & 2 ^ c(-|,«)C (|,« + V&) ^ j. (44) 
The finite temperature case 

As shown in Appendix B, the grand potential at finite temperature can be 
obtained from its expression at T = 0, ljo(s f ), as: 

/OO g x 
o; (/J + Tx) n , 2 (45) 

where the temperature T and the chemical potential fi are given in units of 
the electron mass. If ujq(ij, + Tx) is a slowly varying function of x over a range 
Ax ~ 1, then the integral in Eq. ( f45| ) can be easily expanded in powers of T: 

= u>M + + |^V> 4 V) + • ■ ■ ■ (46) 



In our case of interest, this approximation is valid for the non-oscillatory terms 
4™°', but not for the oscillatory terms, which vary considerably over the range 
Ax ~ 1 (unless the unlikely condition \xT jb <C 1 is met). 

The treatment of the oscillatory terms, which are of the form C 3 (z 1 {(/i 2 — l)/b}), 
for z a negative semi-integer, can be done as follows. Consider placing the 
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leading oscillatory term, 



ui, 



(osc) /// j 2 6 5/2 ^ , ., 
o 



^ 5 VfH. «/»})+- • <«) 



into the integral of Eq. (|4J) . Accordingly, we must evaluate this function at e' F — 
p+Tx andpp 2 = e' F 2 — 1 = p 2 +2fiTx+T 2 x 2 (where T is the temperature in units 
of the electron mass) and weight it with the hump function h(x) = e x / (e x + 1) 2 . 
In the degenerate regime, T <C /i, so the term 0(T 2 ) in the expression for p' F 2 
can be safely neglected. Since the prefactor of the zeta function in fl47| ) is a 
slowly varying function of a;, we can approximate 



m 4 /) 5 / 2 

W (° SC )(T, M ,6)^- — / C (-§, {pl/b + 2x pT/b}) h(x) dx. (48) 

The number of oscillations that fall under the hump will clearly be proportional 
to the factor 2uT/b. We expect the amplitude of the oscillatory magnetization 
to be more or less constant for \xT '/b < 1 and rapidly decreasing for uT/b > 1. 
To obtain an explicit expression we use Hurwitz's Fourier expansion of C(z, q) 
shown in Eq. (|7^) (valid for z < and < q < 1) and integrate term by term 
using the formula: 

e iax h(x)dx= . ™ > ■ (49) 
, smh(7ra) 



Defining the integral in Eq. (|48|) in a generic form as: 

/oo 
C (z, {a + f3x}) h(x) dx, (50) 
-OO 

we find it to be 

/,(«, P ) = (2n)^T(l - z )f,£ + Z7r/2) . (51) 

sinn [2n^pn) 

n—l v 7 

Now, using this result for z = —3/2, the oscillatory part of the grand potential 
becomes: 

::^(T u h\ ~ I m * h*' 2 T V cos^Tmpg/b - tt/4) 1 

W (T, M ,6)_+— 16 _^____^_ 5 _^ + ...|. (52) 

5 The Magnetization 

The T = limit 
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The magnetization (per unit volume) of the electron gas at T = is given by 



Mo 



daj (£ F ,B) 
dB 



2e daj (e F ,b) 

~7 X 



e m 2 d 



2tt 2 db 



db 

&5/2 r i/b ^^M dq 



(53) 



Taking the derivative we find the formal expression for the magnetization 
M (e Fl b) = ?^l b ^ £ F/b H_ 1/2{q) l 2 1 



(1 + bq) 1 / 2 2(1 + bq) 3 / 2 



& 1/2 ^W- 1/2 (p|/&) 



dq 
(54) 



Just as wo was separated into two terms according to the monotonic and the 
oscillating parts of K_i/2, so can we separate the magnetization as Mq(e f , b) = 

M { mon \e Frl b)+M ° sc \e F ,b). The small b expansion for Mo can be obtained 
from Eq. ( |54|) , or directly taking the derivative of the expansion for uuq given in 
Eq. ©: 



M (e F ,b) 



2tt 2 



| - & 1/2 ^ C(-h {Pl/b} + I b cosh- 1 ^)) 



6 3 / 2 (4e 



3 

0(b 3 ). 
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(55) 



The magnetization oscillates as a function of 6 (see Fig. 3), just like the well 
known de Haas-van Alphen effect of non-relativistic electrons in metals. In 
this expression, the oscillation appears in terms of Hurwitz functions of the 
fractional part of p F /b. Notice that, for 6^1, the oscillatory part has an 
amplitude considerably larger than the monotonic part, and even larger the 
larger e F is (see envelope curves in Fig. 3). As such, it could be possible to have 
spontaneous magnetization for sufficiently dense systems at low temperature. A 
discussion about the thermodynamic stability of such state was given in Ref . Q . 
However, at temperatures above some threshold, the oscillation amplitude dies 
out, and with it the possibility of spontaneous magnetization J^.|^]. 

The finite temperature case 

We will be concerned in particular with the physically relevant case where 



Pi/b > 1, 



many Landau levels occupied. In this case, the dominant 



contribution at T = is the leading oscillatory term followed by the leading 
monotonic term, i.e. the second and first terms in Eq. ([sBI), respectively. 
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At finite temperature, the monotonic term becomes: 



M (mon) (T,^,b) 




(56) 



where p 2 = ji 2 — 1. This contribution is small compared to the oscillatory 
part for temperatures smaller than the Landau level splittings, i.e. T < b/fi. 
However, for higher T the amplitude of the oscillations gets smoothed out and 
it is the monotonic part shown above what dominates the magnetization. This 
behavior is shown in Fig. 4. 

On the other hand, the leading oscillatory term at finite T is obtained following 
the same procedure as in the previous section, or directly taking the derivative 
of the grand potential with respect to b: 

Mb-HT,* b) = £ P^lffrtl + • • • • (57) 

k \2b V n smh(47r 7i fi 1 Jo) 

This expression coincides with the results in Ref. p|. We must also remark the 
strong similitude with Landau's result for the non-relativistic case jfjl], which 
is expressed in our notation as: 

MT'CT.M) = - — ^ E ^ 2m ff 2 ~lit (58) 
™' r - v ' 7T ^/2b *r? x Vn smh(47T 2 nT/6) v ' 

This is precisely the non-relativistic limit of our leading oscillating term, shown 



in Eq. ( |57[) : basically, the relativistic chemical potential reduces to fj, = 1 + 
fJ-n.r. — > 1, in units of electron mass, while p F = \J [i 2 — 1 — > y / 2fi n r , denotes a 
non-trivial quantity in all regimes. 

A final important point is the possibility of ferromagnetic behavior, that is, a 
magnetization that is sustained with the magnetic field generated by the same 
system. For a given geometry, the field and magnetization satisfy a relation of 
the form M = jB, where in the case of e.g. a sphere, 7 = 3/87T. A self-consistent 
solution is obtained from this relation and the thermodynamic relation M = 
M(T, fj,, B). This problem has been studied in the past, and our conclusions 
agree with those results. For the self-consistent solution to exist, the oscillations 
are necessary, because the monotonic part of M vs. b is simply too small to 
reach a solution with M — jB. Moreover, for the oscillations not to be thermally 
damped out, the temperature must be below some threshold value ||. Finally, 
there is the question of thermodynamic stability of the self-consistent solution, 
which is at most metastable according to O'Connell and Roussel M. 



Conclusions 
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We have developed a closed analytical approach to solve the thermodynamics 
of a free gas of electrons immersed in an uniform magnetic field of arbitrary 
magnitude. The method is completely relativistic and particularly useful in the 
case of a degenerate gas, a likely situation to be met in the interior of white 
dwarfs, neutron stars, and magnetars. A central role is played by one of the 
least known of the special functions of mathematical physics, the Hurwitz zeta 
function. Its appearance comes about due to a sum over energy levels of the 
form (a + (in) z (Landau levels for the system treated in this paper), labeled 
by a non-negative integer n. As such, the method should be applicable to 
study the thermodynamics of other systems. The grand potential (and therefore 
all thermodynamic quantities) can be expressed as a one-dimensional definite 
integral, which can be explicitly evaluated in several interesting regimes due to 
the various analytic properties of the Hurwitz zeta function. Hence, our work 
provides a unified derivation of several of the results found scattered in the 
literature. We reproduce the de Haas- van Alphen behavior of the magnetization 
in the relativistic gas and the dilution of it at high temperatures. 

We have not included quantum electrodynamics corrections in this treatment. 
Work to extend our results in this direction, in a fully relativistic fashion, is in 
progress. 
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Appendix A: The Hurwitz Zeta Function 

We collect here some of the properties of the Hurwitz zeta function and present 
a (presumably original) derivation of its asymptotic behavior for large q. For a 
detailed account consult, for instance, references |L6| or [ |l7j . 

The Hurwitz zeta function £(z, q) is defined as the analytic extension to the 
whole complex z plane of the series 

cm = m 

71=0 ^ HJ 

valid for Rez > 1 and q ^ 0,-1,-2,.... The resulting function is analytic 
everywhere except at z = 1, where it has a simple pole with unit residue. 



For q ^ one has 



((z,q+l)=t(z,q)-^, (60) 



which, iterated N times, leads to 

N 



E(^j7 - az,q)-az,q + N + l). (61) 

This hnite sum is what appears in the grand potential as a sum over Landau 
levels, and its expression in terms of Hurwitz Zeta functions is what allows us 
to expand the thermodynamic quantities in different limiting scenarios. 

The derivative of C(z, q) with respect to q is again a Hurwitz zeta function: 

£-C{z,q) = -zC{z + l,q). (62) 



Hermite's integral representation 



1 +1 t l^_ z t „_ z+1 f°° sin{z tan- 1 t)dt 



^ - z—^ Z++ 2^ + 2 ^ + J Q (l + ^/ 2(e ^-D ' (63) 

valid for all z ^ 1 and q > 0, can be used to study the large q behavior of £(z, q). 
For large q, the leading contribution to the integral in Eq. (^3|) comes from the 
region of small t. To isolate this contribution we split the integration range 
into [0, a] and [a, oo), where a is a fixed number less than one (for instance, 
a = 1/2). In order to approximate the first integral we shall use the remarkable 
series expansion: 

(1 + tV 2/2 sin(z tan-H*)) = £ ^^k^ > ( 64 ) 



k=0 
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which converges uniformly for \t\ < 1. In ( p4| ) (z)„ is the Pochhammer symbol, 
or shifted factorial, defined by 

(z) n = ^±^ = z(z + l)---(z + n-l). (65) 



Thus, for a given z and < t < a < 1 we have the uniform approximation 

< e, (66) 



sin(ztan 1 t) ^ fc (z) 2 fc+i ^>fc 



i (1 + ^/2 > (2fc+1) , 

where e can be an arbitrarily small positive number and JV = iV(e, z) is a 
sufficiently big number. So we can write 

■ ¥^^v(-if#^ r /7 1 dt+fl( e)> (67) 

o (1 + t 2 ) z / 2 (e 2 ^ t - 1) ^ v ; (2fc + l)!i e 27r «*-l w ' v 7 

with 

l-R(e)l< / ^r4 7 dt<e -— dt = —^. (68) 

Also, up to a correction that vanishes exponentially as q — > oo, we can approx- 
imate 

a .2fc+l poo j.2fe+l n 

e 2nqt _ ! at - J Q e 2n q t_i al \ L > 4 (fe + ±) q 2k+2 ' ^ 

where we have used formula 3.411.2. of reference p8[ . The £>„ in (|6^) arc the 
Bernoulli numbers of even index, i? = lj -B2 = 1/6, -B4 = —1/30, etc. Since 
the integral over [a, 00) in ( |63| ) also gives a exponentially small contribution for 
large q, we finally have the result 

C(z a) = —q- z+1 + -q- z + Q-*- 1 V (-D fc ( z ) 2fc+lg2fc + 2 
U ' 9J z-1 9 + 2 9 +9 ^0 (2k + 2)\q^ 

+0(e)+0(e- aq ). (70) 

Letting e — > (and consequently N — > 00) we obtain the following expansion 
for ((z,q): 

where ( f72| ) follows because the only non-vanishing Bernoulli number of odd 
index is B\ = —1/2. 
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The cases z = — 1/2 and z = 1/2 are of special relevance for the work of this 
paper and we give the corresponding asymptotic expansions explicitly: 

«-§,?) = -^ 3/2 + ^ 1/2 -^- 1/2 + 0(,- 5 / 2 ), (73) 
C(l,q) = -2q^ + 1,-1/2 + l_ q -s/2 + ( q -VZ). (74) 

For general z, the expansion (|7l]) is only asymptotic. It is easy to check that 
its radius of convergence, as a series in the variable 1/q, is zero. However, 
as an aside we should mention that for z = — m, where m is a non-negative 
integer, the series ( |7l| ) terminates, in view of (— m) n — for n > m. Hence the 
series becomes a finite polynomial in q, actually a Bernoulli polynomial, up to 
a multiplicative constant: 

C(-m,q) = 1 B m+1 {q), m = 0,1,2,.-.. (75) 

m + 1 

The Hurwitz zeta function for z < admits the following Fourier expansion in 
the range < q < 1: 

2r(l — z) I /7rz\ cos(27ron) /7rz\ ^ sin(27rgn) \ 

^ q) = x sin (t) s + cos (t) s 4^ • 

V 7 \ 71=1 71=1 / 

(76) 



Appendix B: Finite Temperature from T = 

In this appendix we show that the full finite temperature grand potential shown 
in Eq. (kII) can be obtained from its zero temperature limit as 



il(T,fi) = J n (n + Tx)h(x)dx. (77) 

where h{x) is the Fermi-Dirac hump, 

h(x) = , ^ xo = \ . (78) 

V 7 (e* + l) 2 4cosh 2 (a;/2) V ' 

Actually, result §j% is a special case of the fact that any physical quantity of 
the form 

Q(T, M ) = ^ \{E) I ■ dE (79) 
can be computed in terms of its zero temperature limit Qo(m) as 

Q(T,n)= / Q (n + Tx)h(x)dx. (80) 
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To prove (80) we express the function q(E) in Eq. ( |79[ ) as the derivative of a 
function F(E) defined as: 

F{E)= j E q(E')dE', (81) 
and integrate by parts. After the change of variable E = \i + Tx we obtain 



Q(T,ijl) = J^F(n + Tx) {eX + 1)2 dx. (82) 

But from ((79|) it is seen that the function F(fi) is precisely the value of Q(T, /j,) 
in the limit T — ► 0: 

Qo(/x) = lim Q(T, /i) - /" g(£)dE = F(/i). (83) 



This establishes (|30|). Now, to prove Eq. ( |77| ) for the grand potential, we merely 
notice that the standard expression for it shown in Eq. ( ^l| ) can be turned into 
the form ( |79| ) by expressing the density of states g(E) as the derivative of the 
function G(E) -which is the total number of single particle states up to energy 
E: 

G(E)= j g{E')dE', (84) 

and then integrating by parts to obtain 

n(T, M ) = - jT G(E) ^ E _^ T — ^ dE. (85) 

This expression for the grand potential is precisely of the form shown in Eq. ((79|) , 
with q(E) = —G(E). 

In the degenerate regime, T <C [i — m, the lower limit of integration in ( f77| ) 
can be replaced by — oo with negligible error. Additionally, if the function 
Qo(/U + Tx) varies slowly under the hump (this may not be the case for the 
oscillatory terms; see section ||), then it can be expanded in a Taylor series 
around fj, and the resulting terms integrated one by one. Using now the results 

f°° f°° 7T 2 f°° 77T 4 

/ h(x) dx = 1, / x 2 h(x) dx — — , / x 4 h(x) dx = — — , 

J —CO J —OO J — OG ^ 

/OO 
x 2n h(x) dx = TT 2n | (2 2n - 2) B 2n I , (86) 
-OO 

we obtain the small temperature expansion, 

a(T, = n ( M ) + t 2 q"^) + ^ t 4 ( M ) . . . (87) 
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Appendix C: Reduction of the Oscillatory Con- 
tribution 

Let p(q) be a function in the unit interval [0, 1] periodically extended over the 
real axis, and f(q) an arbitrary function. Then, the integral of the product 
of these functions over an arbitrary interval can be separated into a sum of 
integrals over integer intervals, plus a residual integral, as follows: 



rQ /-i LQJ-i „{Q} 

/ P(q)f(q)dq= / p(q) £ f(q + k) dq + / p(q)f(q+[Q\)dq. 
Jo Jo k=0 Jo 



We now specialize this result for p(q) = ((z, {q}) as the periodic function, 
f(q) = (1 + bq)~ s and Q = p%/b, and use Eq. ([H]) to express the sum in terms 
of Hurwitz functions: 

\j>l/b\-i l [p%/b\-i 

X! fi 4- h(n -I- h\] S ~ h» 



,_ u [l + b(q + k)] s b* ^ [(q+l/b) + k] s 

= I [C («. 1 + V&) - C (s, g + 1/6 + Lpp/H )] • 
We thus find the integral expression 

= ,r' /j , (88) 



o (! + &<?)* * 6 s W (q + l/6+b?/6J) 



+ ^ ^Q)[C(s,q+l/b)-C(s,q+l/b+[pl/b\)] dqj. 

Using this generic result one derives the expression given in Eq. (|3^) for the 
oscillatory piece of the grand potential. 

Appendix D: Small b Expansions 

In this appendix we show that 



pl/b rj-l a +i) 9 1 



-a- 5 -)- — 

L5 U 2> 126 
+ 0(6 3/2 )... (89) 



and 

r cJ 7mr dq - -§<(-»- i5«-i>»-5<H)" 

+ |fcH,04/t}) + ^^c(-f,04/t}) 



2 



+ |^C(-|,{^})+0(6 3 )... (90) 
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The first expansion is needed in order to find the small b behavior of the non- 
oscillatory piece of the grand potential, Eq. (34), and the second expansion is 
required for the oscillatory piece, which is exhibited in Eq. ( |36| ) already making 
use of the formula (pq). 



To prove Eq. 



we first define the integral 
R(b) = 



' PF/b C 3 H,g + i) 



VI + bq 



dq, 



and then try to expand it in powers of b. The first two terms of the expansion 
are easily found, using the results: 



flWUo = 5«-5) + 60 . 



6=0 



2 
3 
2 
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1 



^C(-I) 



1260' 



which follow from the facts that <^3 ( — ^ , 9 + 1) decreases like q 5 / 2 when q — > oo, 
and that the antiderivative of Ca(— \, q), which is |^3(— |, g), vanishes like q~ 3 / 2 
when q — > oo . Indeed, we can directly calculate at b = 0: 



R(b)\ 



6=0 



(3(~h,q)dq 



2 - 
3 



CaC-f, 1) 



-C(- a ) + - 

3 U 2 > 60' 



(91) 



and also 



6=0 



lim 

6^0 



l+pl/b) 



£ F b 2 



Ca(-i« + l) 



(l + bq)V- 



dq 



l)dq. 



(92) 



In this last expression, the 6—^0 limit of the first term inside the square bracket 
vanishes due to the large-q behavior of Ca(— |> ?)) while the last integral can be 
done by parts: 



1 

2 

1 2 
3 ' 5 



Cs(-1,9 



l)dq 



Cs(-|,i) = ^c(-l)- 



1260 



We should point out that the next term in the expansion J89| ) is not 0(b 2 ), as 
it would be in a regular Taylor series, because R"(0) does not exist. To get this 
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singular term one should go back to the original definition of R(b) and subtract 
the leading term in the asymptotic expansion of Cs{—^,q)'- 



i 



1920 g 5 / 2 



+ &(-!, 9), 



(93) 



where Ci(~ |j l) is °f order 0(q 9 / 2 ). Then, the resulting term in the integrand, 
(1 + g)~ 5 / 2 (l + fog)" 1 / 2 , can be integrated explicitly, and its small b expansion 
reads 



Pl/b 



1 



„ , T -=dg=|-|&+|5(2^-l) b*' 2 + 0{b 2 ) (94) 
o (1 + g) 5 / 2 v! + h 3 3 3p£ 



Therefore, the 0(& 3 / 2 ) term in the expansion of R(b) is actually 

L_ £l C2» 2 - 1) b 3/2 

2880 p? 1 Pf j b ■ 

With this result we complete the proof of Eq. (|S9]). 



(95) 



We now proceed to prove Eq. (^0|), which is relevant for the small b expansion 
of the oscillatory piece of the grand potential. First we must realize that the 
integral in Eq. (|o|) is of the type studied in Appendix C, and in particular the 
result shown in Eq. ( |S8| ) applies. We will then proceed to expand that result 
for small b. To expand the first integral in Eq. (pq), we just need to use the 

binomial expansion of (g + ^ + \p 2 /b\ ) , where q <C 1/6 + [p 2 /b\, and then 
use the integrals: 



Q a-hi)dq = |cc-i,o) - |cc-i), 



(96) 



J qa-hl)dq = -^C(-f,Q) + ^C(-§) + ~QC(-§,Q). (97) 
We thus find the expansion for the first integral of Eq. (^8|) : 



( q + i/b+[p 2 /b\) 



1 



{i + b[ P i/b\y/ 2 
b 

' (l + b[ P 2 Jb\)V 2 



|c(-|,{rf/6})-|c(-|) 



+ 0(b 2 )... (98) 



Notice that the terms of the form (1 + b \pl/b\) n l 2 £(— TJ -^) will cancel in 
the full oscillatory piece, Eq. (p6|), when expansions (p8|) and (102) (see below) 
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are combined. For the remaining terms containing (1 + 6 Ipf/M) in the de- 
nominator we write 1 + 6 [Pf/^J = (1 +-Pf) — b{p F /b} and perform a binomial 
expansion, which leads to a further cancellation of all the terms of the type 
{p 2 F /b} k ((z, {pl/b}) with k > 1. 



Now we need to expand the second integral in Eq. (|88|). In this case we must 
expand the £ functions inside the squared brackets for large values of their 
argument. We thus use the large-g expansions of £(±i,g) shown in Eqs. ( |73|) 
and (|74|), together with the definite integrals [|l9| 

l 

C(-i«)dg = 0, (99) 
2 

3 C 



qC{-\,q)dq = -£(-§), (100) 



q 2 C(-hq)dq = |C(-|) - ^C(-|). (101) 



We thus get the expansion: 

i r 1 



VbJo 



C(-|,g) [C(|,g + - C(|, q + l/b+ [p%/b\)} dq 



+ 0(6 2 ).. 



1 - 



(1 + &LpI/&J) 1/2 



1 - 



(1 + &LPf/&J) 3/2 . 

(102) 



which completes the proof of Eq. ( |90| ) up to order b. The 0(6 2 ) terms are 
obtained in a similar fashion. 

Finally, and for completeness, we also present here the small b expansion of the 
explicit terms in (|34|), which reads 

l(l-6+l6 2 )lnf £F + ^ + 
2 V 6 V 1 + ^ 



- P l) 1/2 e F -b^ 



(b + plf 2 



e F - 6 3 / 2 

1 

2 



-6 [cosh 1 (e F )-p F e F ] 

1 u-l/ x 1 1 - 1 

= -cosh (e F ) - -p F e F + -? 



3^F e F + — 6 2 cosh 1 (e F ) 

60 1260 1 ' y ' 

One should notice that all the terms that are independent of the Fermi energy 
(i.e. p F or e F ) are spurious, cancelling between the different expansions and 
thus leading to result (|38|). 
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Figure 1: The particle density at T=0 as a function of the Fermi energy, for 
fixed magnetic field b = 0.1 (smooth curve) and b = 0.5 (bumpy curve). 
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Figure 2: The particle density at T=0 as a function of the magnetic field, 
fixed Fermi energy e F = 1.3 (lower curve) and e F = 1.5 (upper curve). 
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Figure 3: The magnetization at T 
magnetic field b, for Fermi energy e F 
envelopes of the curve, and the upper 
and 16). 



which oscillates as a function of the 
4. Also shown are the upper and lower 
velopes for larger Fermi energies (e F = 8 
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Figure 4: The upper envelope of the magnetization as a function of the temper- 
ature, for fixed chemical potential fi = 100, and magnetic field b = 0.1, 1 and 
10 (lower, medium and upper curve, respectively). 
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Figure 5: The upper envelope of the magnetization as a function of the magnetic 
field b, for two temperatures T — 10~ 5 and 1CP 7 (lower and upper curve, respec- 
tively). The chemical potential is fixed at /i = 100. The dashed line represents 
the macroscopic relation M = -S-B, for a uniformly magnetized sphere. 
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Figure 6: The functions ((z, {q}), for z = —1/2, —3/2, —5/2 (large, medium and 
small amplitude, respectively). 
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